Scale decomposition of molecular beam epitaxy 

Z. MoktadiiE 

In this work, a study of epitaxial growth was carried out by means of wavelets formalism. We 
showed the existence of a dynamic scaling form in wavelet discriminated linear MBE equation where 
diffusion and noise are the dominant effects. We determined simple and exact scaling functions in- 
volving the scale of the wavelets when the system size is set to infinity. Exponents were determined 
for both, correlated and uncorrelated noise. The wavelet methodology was applied to a computer 
model simulating the linear epitaxial growth; the results showed a very good agreement with an- 
alytical formulation. We also considered epitaxial growth with the additional Ehrlich— Schwoebel 
effect. We characterized the coarsening of mounds formed on the surface during the nonlinear phase 
using the wavelet power spectrum. The latter have an advantage over other methods in the sense 
that one can track the coarsening in both frequency (or scale) space and real space simultaneously. 
We showed that the averaged wavelet power spectrum (also called scalegram) over all the positions 
on the surface profile, identified the existence of a dominant scale a*, which increases with time 
following a power law relation of the form a* ~ i", where n ~ 1/3. 

PACS numbers: 



I. INTRODUCTION 



Over the last two decades, many aspects of molecu- 
lar beam epitaxy (MBE) were theoretically investigated 
within a phenomenological framework [U, S B 0| • Phe- 
nomenological continuum models consider the surface of 
the growing film as a continuous variable of the posi- 
tion where overhangs are not allowed. They have the 
credit to explain many aspect of the surface morphol- 
ogy of growing films 0, The MBE process can be 
described as follow: atoms are deposited on the film sur- 
face from the gas phase, where they undergo a thermally 
activated diffusion or desorption back to the gas phase. 
Once absorbed, atoms can combine to form a dimmer, or 
attach to the steps of existing islands on the surface. A 
whole atomic layer is completed once all islands on the 
surface have coalesced. Smooth surfaces are obtained 
in a layer-by-layer growth mode in which a new layer 
starts to form only when the layer underneath is fully 
grown. However, experiments provide evidence that the 
layer-by-layer growth mode do not occur in many situa- 
tions (see for example references Q, 0] and Q). This 
ideal situation is suppressed by two dominant effects: 
shot noise and instabilities that arise from the so-called 
Ehrlich— Schwoebel (ES) effectQ. Shot noise originates 
from different mechanisms such as deposition, diffusion 
or nuclcation. The ES effect is due to the asymmetry in 
attachment— detachment kinetics across an atomic step, 
where atoms have to overcome an energy barrier when 
descending the step. This triggers an ascending atomic 
current, which is responsible for a morphological instabil- 
ity. In this case, the amplitude of a small perturbations 
on the flat surface will exponentially increase in time. 
This instability can be balanced by the introduction of a 



stabilizing mechanism such as MuUins-like current arising 
from thermodynamic relaxation through surface diffusion 
[iQ] or from fluctuations in the nuclcation process of new 
forming islands [ll|, • The ES effect and diffusion cur- 
rents will induce the formation of a mound-like structure 
on the surface which coarsen as time progresses. 
A phenomenological continuum model describing the sur- 
face growth incorporating the two conserving mecha- 
nisms mentioned above, can be formulated in one dimen- 
sion as: 



dh 
'dt 



(1) 



where h is the surface height, jd is the ES destabiliz- 
ing current, jg is MuUins stabilizing current and ri{x,t) 
is noise function representing the stochastic character of 
the growth. This function is assumed to be a Gaussian 
white noise with zero mean, or a spatially correlated noise 
(long-range correlations), i.e.: 

< r;(x,t)77(x',<') > = 2F(5(x-x')5(t-t') 

< r;(x,t)77(x',<') > - 2F|x-x'|2p-i(5(<-<') (2) 

where F is a constant and p (0 < p < 1/2) is an expo- 
nent characterizing the decay of spatial correlations. The 
Fourier transform of the noise correlators above is given 
by: 

<77(q,t)?7(q',t') > = 2F(5(q + q')<5(i - i') 
<77(q,<)77(q',t') > = 2Dpq-^PS{<i + c^)S{t - t') (3) 

The prefactor Dp is given by: 

D..£r.-'c„w...£Q^ (4) 
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Where F is the Gamma function. 

A simple model for the currents jd and js can be ex- 
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pressed as [11] : 



A. The linear MBE equation in wavelet domain 



jdim) = 1 - 



Too 



(5) 



where TO = §f is the surface slope. The parameters v 
and K are positive constants related to microscopic pro- 
cesses of deposited atoms on the surface [13|. The pa- 
rameter Too is the so-called "the magic slope". The form 
of the destabilizing current jd predicts slope selection: 
i.e. the mound slopes will asymptotically reach a con- 
stant value Too. Equation ^ has been investigated by 
mapping it to the phase ordering problem or by map- 
ping it to a one dimensional system of interacting kinks 
p^ . The scenario predicted by ^ is the following: The 
competition between ES effect and the surface diffusion 
will lead to a mound like structured surface with a well 
defined period. Later in time, the mounds will coarsen 
because of the nonlinearity of the current jd- The period 
of the mounds X{t) scales with time as X{t) ~ t" where 
n= 1/3 

The purpose of this a paper is to characterize the MBE 
process described by ([T]) using wavelets formalism. Two 
cases are considered here: linear stable growth where only 
surface diffusion and noise are taken into account, and 
growth incorporating the additional effect of ES barrier. 
In the former case, the scaling functions and exponents 
are derived in wavelet space in the presence of both cor- 
related and uncorrelated noise. In the latter case, the 
coarsening is discriminated through wavelet decomposi- 
tion of the evolving surface profile and characterized by 
the wavelet power spectrum. The advantage of using 
wavelets is that one can track the coarsening process in 
real space and the frequency space simultaneously. 
This paper is organized as follow: we first consider the 
linear stable MBE growth(section II) by performing an 
analysis of growth in wavelet space. Section III is de- 
voted to nonlinear unstable MBE process. We close the 
manuscript with a conclusion. 



II. THE LINEAR MBE 

In this section we will use the wavelet formalism in 
order to determine the properties of the linear MBE 
through the scale discrimination of the growth process, 
that is, the determination of the scaling function and ex- 
ponents as a function of the scale when the system size is 
infinite. The wavelet transform of a profile h{x) is given 
by [11: 



T (a, x) 



1 



y-x 



a 



dy 



(6) 



where a (a > 0) is the scale parameter and is the 
mother wavelet. The transform T{a, x) is a scale-position 
decomposition which expands a function h{x) in wavelets 
basis, whose elements are constructed from a single 
mother function: the mother wavelet. 



When diffusion is the dominant process in the surface 
dynamics and all other destabilizing effects are neglected, 
(fT|) is reduced to a fourth linear form (for simplicity we 
set K=l): 



dh ^. , , 

= -V^h{x) +r]{x,t) 



(7) 



This equation can be transformed by the help of Hermi- 
tian wavelets [l3, EBl which are the successive deriva- 
tive of the Gaussian function exp(— Hermitian 
wavelets u„(x) (here n is the derivative's order) have re- 
cursive properties allowing the transformation of linear 
equations from direct space to wavelet space. Applying 
the wavelet definition ([6]) to ([7]) we can write: 



d_ 

di 



+S{k, X, t) 



(8) 



Where: 

Sn{K,x,t) ^ yfn, I f]{y,t)uniK iy - x))dy (9) 



is the wavelet transform of the noise ri{y, t) and k ~ ^/ a, 
where a is the scale. Now, Integrating by parts the in- 
tegral ([8]), and use the recursive properties of Hermitian 
wavelets (see reference [l^), we arrive at: 

d ( d"^ d \ 

-Q-^Un = -5l(«) ( ^^""^ ~ 5^"" j ^ 92{K)Un + S{K,X,t) 

(10) 



where: 



52(«) -«;'(n+ 1/2)2 



(11) 



Equation (fTO|) describes the linear MBE growth in 
wavelet space. This equation will not be subject of fur- 
ther analysis in the present paper. 



B. Dynamic scaling 

In a previous work ^2^1 we applied wavelet formalism 
to Edwards- Wilkinson equation. This investigation was 
concerned with dynamic scaling in terms of scale a but 
not the system size. We derived an exact and simple ex- 
pression for the scaling function. The dependence of the 
surface width a was found to be a scaling law of the form 
a{a,t) ~ a for uncorrelated noise and a{a,t) ~ a^^^ for 
correlated noise. Here, we will apply the same formalism 
to ([7]). The the lateral system size L is taken to be in- 
finite and therefore the dependence on L is suppressed. 



3 



Since each decomposition T{a, x) does not result in a sine 
wave, it is possible to calculate its power spectrum. By a 
simple change of variable in ([6]), we compute the Fourier 
transform of each decomposition : 

f{a,q,t) = ~V^h{q,t) y^ACOe^'^de 

= -^h{q,t)i^{~aq) (12) 

where h and ip are the Fourier transforms of the height 
h and the mother wavelet ip respectively. The power 
spectrum at a scale a is then : 

r{a,q,t) = <f{a,q,t)f{a,-q,t)> 

= aP{q,t)\^{~aq) p (13) 

The function ipiq) is the Fourier transform of the mother 
wavelet. For the second order Hermitian wavelet ^/'(a;), 
the function i/;(g) is given by: = \/27rq^ exp(— 9^/2). 
The quantity P{q) is the power spectrum of the surface 
height. By Fourier-transforming ([7]), we obtain the solu- 
tion for a flat initial condition: 

h(q,t) = / e'i'^'-'^{q,T)dT (14) 
Jo 

where h {q, t) and i) (q, t) are the Fourier transforms of 
the height h, and the noise rj respectively. 
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FIG. 1: The scaling function f{x). 
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TABLE L Values of the exponents in the two cases: the EW 
equation and the linear MBE equation in one dimension. Be- 
tween brackets are values of the exponents observed in the 
scaling with the system size. 



1. Uncorrelated noise 



Using (I14|) and ([3]), the power spectrum is given by 
< h{q,t)h{—q,t) >, i.e.: 



P{q) = 2F 



(1 - exp(-2g4t)) 



(15) 



We can now calculate the surface width at each scale. 
From we get: 



(7{a,t) — / T{a,q,t)dq 
Jo 

= a^Fn [ 
Jq 



l-e-"« M e-" « dq (16) 



The upper cutoff qmax is of the order of the inverse lat- 
tice constant; we assume that the correlation length is 
larger than 1/qmax and we set qmax to infinity in 
Performing the integration we arrive at: 




order n. Thus, the dynamic exponent is z = 4. The satu- 
ration value of the width scales as a sat oc a", with a = 2. 
The scaling function f{x) is not defined at a; = 0, but 
has the asymptotic limit f{x) oc x^^"^ for a; << 1 ( this is 
easily shown by performing a series development to the 
second order of f{x)) and f{x) \/2'k^/'^ for x » 1. 
This function is displayed in figure [TJ In one dimen- 
sion, the scaling law dependence of the surface width at 
a scale a, involves integer exponents in the two cases, the 
EW [l3] and the linear MBE equations, unlike the scaling 
exponents observed in the dependence with the system 
size, which are fractional. Table [T] summaries the value 
of exponents in both cases in one dimensional space. 



2. Correlated noise 

Similar to the above calculations, we will determine 
a different set of exponents by computing the surface 
width corresponding to each wavelet scale a taking the 
substrate size infinite, in the case of spatially correlated 
noise. In this case the power spectrum is given by: 



Pp{q) = 2Dpq 



(17) 



We have for the surface width at the scale 



where f{x) = ^2^^ - ^ exp(^)if (1/4, l/(16x)) 
and K{n, x) is the bessel function of the second kind of 



{a,t) 



= 2a''Dp7r 



(18) 

2 2 

" dq 
(19) 
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In this integration we have taken the upper cutoff to be 
infinite. By simple change of variable we have: 



(Tp{a,tf = TrDpa^+^P 



-p-l/2 



(1 



here, 



P \ ,A 



(20) 



is a scaling 
^/x for cc << 1 and 



function which has the limit fp{x) 
fp{x) « A/r(l/2 — p), for x >> 1, where T is the gamma 
function. The result for the uncorrelated case is retrieved 
for p = 0. Thus, for correlated noise, the roughness ex- 
ponent is a(p) = 2 + p while z = 4, independent of p. 
For the EW growth model, the roughness exponent was 
found to be uew = 1 + p [2(| |. 
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C. application to a computer model of the linear 
MBE 

In this section we will apply the results obtained above 
to a computer model that simulate the linear molecular 
beam epitaxy in the case of uncorrelated noise. This 
model was developed in [2l| to simulate the linear MBE 
growth. In this model, atoms are randomly deposited on 
a linear substrate and undergo a diffusion to neighboring 
sites in order to maximize their curvature k. If the height 
at a site i is h{i), then the deposited atom at this site 
will move to a site j with the maximum value oi k = 
hj^i — 2hj + hj-i. The diffusion length I is such that 
i — l<j<i + l. Simulation is carried out for I = 2 and 
a substrate size of L = 10^ sites. In figure [2] we show 
the plot of the surface width calculated at wavelet scales 
a =2,3 and 4 by performing the wavelet transform of the 
simulated profile and by using the expression a{a,t) — 
{{T{a, i) - T(a))2)^ where T(a) is the spatial average of 
T{a,i), i = 1..N. As can be seen, there is an excellent 
agreement between the simulated value of cr{a, t) and the 
scaling form (fT7|) . The inset of figure [5] shows a good 
data collapse confirming the scaling form (fT7|) . 

III. GROWTH WITH EHRLICH SCHWOEBEL 
BARRIER 

The current form j„i in © represents a continuous 
uphill current leading to a slope selection i.e the mounds 
slope converges to a constant value toq. This current is 
countered by the diffusion current jd preventing an in- 
definite increase of the mounds height. A simple linear 
stability analysis of ([1]) shows that this growth process is 
unstable against small perturbations with wavenumbers 
smaller than a critical value kc = \J v j K. This implies 
that the initial growth stage is characterized by the for- 
mation of mounds on the surface with the typical size 
A — 27r-\/2/gc- After this initial phase, non-linearities be- 
come relevant triggering the coarsening of mounds. The 



FIG. 2: Computer simulation results showing the evolution of 
the surface width for three different scales a=2 (•), a=3 (■) 
and a=4 (a). Full lines are the fit to the analytical expression 
(|17p . The inset shows the data collapse to a single curve 
confirming dynamic scaling form (|17|l . 



scaling hypothesis implies that the coarsening behaviour 
is statistically self-similar i.e. surface mounds are similar 
in time domain up to a scaling by the average mound size 
A. Under the scaling hypothesis, structure functions such 
as the height-height correlation function can be writing 
as C[r) — a{t)^6{r / \{t)), where a is the surface width 
and is a scaling function. The evolution of A and a is 
expressed as: 

A - 

a ^ (21) 

The current form given by ([5]) leads to slope selection 
and the value of the scaling exponent n = | [1]. 
We will show in the following, that the coarsening pro- 
cess can be well characterized by wavelets formalism. 
The advantage of this formalism is that one can track 
the coarsening in both, scale (or frequency) and the spa- 
tial position simultaneously. In addition we can quantify 
the coarsening we determine the evolution of a quantity 
called the scalegram 17] which is the counterpart of the 
power spectrum in Fourier analysis. 

We first start by solving ^ numerically. The wavelet 
transform of the generated profiles is then computed at 
different times. The wavelet's power spectrum is defined 
as the squared magnitude of the wavelet transform which 
is the analog of the power spectrum in Fourier analysis. 
Figure [3] displays the time evolution of the wavelet power 
spectrum at t = 100, 500 and 1000 for v/K = 2 and 
N — 500; the result is averaged over 100 independent 
runs. The Hermitian wavelet of order 1 is used in these 
calculations. We can clearly notice the pattern formed at 
the early stage of growth where the power is concentrated 
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FIG. 3; Color map plots showing the time evolution of the 
wavelet's power spectrum of the evolving profile for three dif- 
ferent times t = 100, 500 and 1000. The early stage {t = 100) 
shows the power concentration in a small band of scales sur- 
rounding the dominant scale as a result of linear instability. 
As time advances, this power is spread through larger scales 
as a result of the coarsening of mounds. The color scale is in 
arbitrary units. 



around the dominant scale as a result of the linear insta- 
bility. As time advances, the power shifts and spreads 
through higher scales, indicating the coarsening process. 
In general, the mounds coarsening is characterized by 
the determination of the lateral correlation length or by 
determining the wavenumber corresponding to the max- 
imum of the Fourier power spectrum. The latter is re- 
lated to the mound size A via the relation qmax = 27r/A. 
One can efficiently characterize the coarsening process 
by performing the calculation of the the scalegrampTf. 
The scalegram is the spatial average of the wavelet power 
spectrum: 



S{a,t) =< \T{a,x,t)\'' >, 



(22) 



The main advantage of the scalegram over the Fourier 
power spectrum is the fact that only few realizations are 
required to obtain an accurate scalegram. We can show 
that the scalegram defines a dominant scale at which it 
reaches its maximum value. We can obtain an analytical 
form of S' as a function of the scale a and the mound size 
A = X{t). Using the wavelet definition ([6]) we have: 

S{a) = [[ <h{x + ^a)h{x + (:a)>^^P{^)'ilj{Od^dC 



(23) 



The kernel C(^,C,a) ee< h{x + £,a)h{x + (a) > is the 
two points correlation function which depends only on 
the difference | — ()a \ for the present process, that is: 



FIG. 4: Plots of the scalegram obtained by numerical simula- 
tion(see text) at three different time t =1000, 1500 and 3000. 
Continuous lines are fits to the analytical form (|25p . 



C(C,C,«) = CIL(?- 

form for C (r) [22] : 



()a I). We will choose the following 



C{r) = exp 



r 



2'Kr 



(24) 



where a is the surface width and I is a parameter char- 
acterizing the decay of the correlations. The lateral cor- 
relation length is defined as C(r) = cr^/e and it is a 
function of both I and A. To obtain a simple analyti- 
cal form of S*, we use the Hermitian wavelet of order 1 
( ■0(x) = — a;exp(— a;^/2). Integrating ([23| in (— oo,-|-oo) 
and using ([24|) . we have: 



S{a) = /(a. A, I) exp 

Where: 

/(a,A,0 



(4a2 +p)A2 
47raV2Z(27r2z4 + /2;^2^4„2^2) 



A2(4a2 + P)5/2 



(25) 



(26) 



We determined the scalegram from numerical simula- 
tions of (U]) by computing the wavelet transform of the 
obtained profile and using the definition (f22|) . Figured] 
shows the results at three different times t =1000, 1500 
and 3000. The Hermitian wavelet of order 1 was used 
with a system size of = 500. The ratio v/K was set 
to 1.5, and the result was averaged over 100 independent 
runs. As can be seen, the agreement between numerical 
solution and the analytical form of the scalegram is ex- 
tremely good. Equation ((25)) predicts the existence of a 
dominant scale a*, which maximizes S, proportional to 
A. Thus, the dominant scale evolves following the same 
power law followed by A i.e. (I^T]) . 
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FIG. 5: The evolution of the dominant scale (corresponding 
to the maximum of the scalegram) obtained from numerical 
simulations, showing a power law scaling with an exponent 
n ~ 1/3. The inset shows the same result in a log-log plot. 



This result is not surprising since the continuous 
wavelet transform gives the information on to what ex- 
tent the frequency content of the analyzed profile, in the 
neighborhood of an arbitrary position a;, is close to the 
frequency content of the wavelet at a given scale. To 
check the validity of the above we numerically computed 
the value of the dominant scale as a function of time. 
This is shown in figure The power law nonlinear regres- 
sion fit is also shown with an exponent n — 0.31 ± 0.01, a 



value consistent with the coarsening exponent n = 1/3 . 
A log-log plot of the result is also displayed at the inset 
of figure [5] confirming the power law scaling. 



IV. CONCLUSION 

In conclusion we investigated the epitaxial growth us- 
ing wavelets formalism. Two cases were considered: the 
linear case where only atomic diffusion is the dominant 
process and the case where both atomic diffusion and 
Ehrlich— Schwoebel barrier are the dominant processes. 
In the former case, the linear equation is decomposed us- 
ing a wavelet filter, allowing the discrimination of growth 
dynamics at each scale. We determined the scaling func- 
tions of the width corresponding to each wavelet decom- 
position of the surface profile, in two cases: growth with 
uncorrelated noise and growth with correlated noise. An- 
alytical results were compared to a computer model simu- 
lating the linear growth in the case of uncorrelated noise. 
A good agreement was found between theory and com- 
puter simulation. 

Growth incorporating ES effect alongside atomic dif- 
fusion is characterized by numerically computing the 
wavelet power spectrum. Coarsening process is quanti- 
fied by the so called the scalegram, which revealed a time 
dependent dominant scale where the scalegram reaches 
its maximum. The dominant scale is proportional to the 
mound size i.e. following a power law with an exponent 
n ~ 1/3. An analytical form of the scalegram is de- 
termined as a function of the scale and the mound size. 
This analytical form was compared to numerical results 
showing a good agreement between the two. 
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